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Abstract. Drained thaw lake basins (DTLB’s) are the domi¬ 
nant land form of the Arctic Coastal Plain in northern Alaska. 
The presence of continuous permafrost prevents drainage and 
so water tables generally remain close to the soil surface, cre¬ 
ating saturated, suboxic soil conditions. However, ice wedge 
polygons produce microtopographic variation in these land¬ 
scapes, with raised areas such as polygon rims creating more 
oxic microenvironments. The peat soils in this ecosystem 
store large amounts of organic carbon which is vulnerable 
to loss as arctic regions continue to rapidly warm, and so 
there is great motivation to understand the controls over mi¬ 
crobial activity in these complex landscapes. Here we report 
the effects of experimental flooding, along with seasonal and 
spatial variation in soil chemistry and microbial activity in 
a DTLB. The flooding treatment generally mirrored the ef¬ 
fects of natural landscape variation in water-table height due 
to microtopography. The flooded portion of the basin had 
lower dissolved oxygen, lower oxidation-reduction potential 
(ORP) and higher pH, as did lower elevation areas through¬ 
out the entire basin. Similarly, soil pore water concentra¬ 
tions of organic carbon and aromatic compounds were higher 
in flooded and low elevation areas. Dissolved ferric iron 
(Fe(III)) concentrations were higher in low elevation areas 
and responded to the flooding treatment in low areas, only. 
The high concentrations of soluble Fe(III) in soil pore wa¬ 
ter were explained by the presence of siderophores, which 
were much more concentrated in low elevation areas. All 
the aforementioned variables were correlated, showing that 
Fe(III) is solubilized in response to anoxic conditions. Dis¬ 
solved carbon dioxide (CO2) and methane (CH4) concentra¬ 
tions were higher in low elevation areas, but showed only 
subtle and/or seasonally dependent effects of flooding. In 
anaerobic laboratory incubations, more CH4 was produced 


by soils from low and flooded areas, whereas anaerobic CO2 
production only responded to flooding in high elevation ar¬ 
eas. Seasonal changes in the oxidation state of solid phase 
Fe minerals showed that net Fe reduction occurred, espe¬ 
cially in topographically low areas. The effects of Fe re¬ 
duction were also seen in the topographic patterns of pH, 
as protons were consumed where this process was prevalent. 
This suite of results can all be attributed to the effect of wa¬ 
ter table on oxygen availability: flooded conditions promote 
anoxia, stimulating dissolution and reduction of Fe(III), and 
to some extent, methanogenesis. However, two lines of ev¬ 
idence indicated the inhibition of methanogenesis by alter¬ 
native e- acceptors such as Fe(III) and humic substances: 
(1) ratios of CC>2:CH4 evolved from anaerobic soil incuba¬ 
tions and dissolved in soil pore water were high; (2) CH4 
concentrations were negatively correlated with the oxidation 
state of the soluble Fe pool in both topographically high and 
low areas. A second set of results could be explained by 
increased soil temperature in the flooding treatment, which 
presumably arose from the increased thermal conductivity of 
the soil surface: higher N mineralization rates and dissolved 
P concentrations were observed in flooded areas. Overall, 
these results could have implications for C and nutrient cy¬ 
cling in high Arctic areas where warming and flooding are 
likely consequences of climate change. 


1 Introduction 

The height of the water table controls the balance of anaer¬ 

obic and aerobic processes in wetlands and other periodi¬ 
cally flooded systems (Blodau, 2002; Reddy and DeLaune, 
2008; Kogel-Knabner et al., 2010). Surface water restricts 
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the diffusion of oxygen (O2), creating anoxic conditions that 
promote anaerobic processes such as methane (CH4) pro¬ 
duction or anaerobic respiration using a variety of alterna¬ 
tive electron (e _ ) acceptors, such as nitrate (NO^), man¬ 
ganese (Mn(IV)), iron (Fe(III)), humic substances and sul¬ 
fate (SC>4 _ ). Surface water can also have other important 
physical effects on soils; for example in soils with thick sur¬ 
face organic layers, increased water content greatly increases 
thermal conductivity (Kujala et al., 2008; O’Donnell et al., 
2009), and so flooding or ponding can lead to deeper heat 
penetration and active layer depth in soils with permafrost 
(Brown, 1967; Shiklomanov et al., 2010). 

Northern permafrost soils contain vast reservoirs of carbon 
(C), and could be a major source for carbon dioxide (CO2) 
and CH4 in the next century as a result of climate warming 
(Schuur et al., 2008; Tarnocai et al., 2009). In the past cen¬ 
tury, high latitude sites have warmed faster than the rest of 
the planet (IPCC, 2007), and this warming trend has already 
caused significant increases in CH4 production from Arctic 
ecosystems (Bloom et al., 2010). However the C balance 
of Arctic ecosystems has shown complex responses to cli¬ 
mate change over the past few decades (Oechel et al., 1993; 
Oechel et al., 2000). Increased plant productivity could help 
offset losses of soil C, particularly if rates of N mineraliza¬ 
tion increase in response to warming (Shaver et al., 1992). 
On the other hand, higher N availability can also stimulate 
decomposition and lead to further losses of soil C (Mack et 
al., 2004). The responses of any given Arctic ecosystem to 
climate change, in particular the rates of ecosystem CH4 and 
CO2 production, will likely depend on local factors, such as 
microtopography and soil properties, and how they influence 
the responses of nutrient cycling rates and the availability of 
alternative e“ acceptors to changes in temperature and hy¬ 
drology. 

The water table in many high Arctic ecosystems is af¬ 
fected by complex microtopographic patterns that arise from 
the formation of ice wedge polygons (Brown, 1967). Over 
many freeze-thaw cycles, ice wedges form in the soil, push¬ 
ing up regions that become the rims of low-centered poly¬ 
gons, causing water to accumulate in the lower polygon cen¬ 
ters. The ice wedges under the rims of these features can 
melt and subside, creating high centered polygons (Billings 
and Peterson, 1980). The presence of this complex spatial 
pattern complicates predictions of the ecosystem response to 
climate change. Previous studies have shown that the net ef¬ 
fect of water- table manipulations on the C cycle depend on 
the plant community and location on the landscape (Upde- 
graff et al., 2001; Strack et al., 2006; Stack and Wadding- 
ton, 2007). Furthermore, climate change will have complex 
effects on the hydrology of northern soils. As northern lati¬ 
tudes have warmed, areas of discontinuous permafrost have 
degraded, leading to draining of lakes; however in areas of 
continuous permafrost, lake coverage has increased as a re¬ 
sult of themokarst activity (Smith et al., 2005). Predicting 
ecosystem responses to climate change in the Arctic requires 
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an understanding of how water table and microtopography 
control biogeochemical cycling and microbial processes such 
as respiration and methanogenesis. 

Drained thaw lake basins (DTLB’s) make up a significant 
portion of the landscape of the Arctic Coastal Plain (Billings 
and Peterson, 1980; Hinkel et al., 2003; Bockheim et al., 
2004). These features originate as lakes that drain due to 
various erosional processes and become vegetated. As these 
DTLB’s develop over the course of hundreds to several thou¬ 
sand years, ice wedge polygons create complex topography 
within these basins. One such DTLB near Barrow, Alaska 
was the location of a large scale water-table manipulation 
(“Biocomplexity”) experiment (Zona et al., 2009; Olivas et 
al., 2010). Previous work at the Biocomplexity site showed 
that much of the basin is prone to anoxia, and that dissim- 
ilatory Fe(III) reduction is a particularly important anaero¬ 
bic process that contributes significantly to the ecosystem 
C budget (Lipson et al., 2010). This process is generally 
thought to competitively inhibit methanogenesis by provid¬ 
ing a more thermodynamically favorable respiratory pathway 
(Roden and Wetzel, 1996; Blodau, 2002; Bond and Lovley, 
2002; Kiisel et al., 2008; Jerman et al., 2009). And while 
CH4 production has been studied in the Arctic Coastal Plain 
near Barrow, AK (Rhew et al., 2007; Zona et al., 2009; von 
Fischer et al., 2010), no studies have yet explored how Fe(III) 
reduction interacts with CH4 production in this ecosystem. 
Additionally, our understanding of Fe(III) reduction in this 
ecosystem has yet to be integrated with a study of fine-scale 
spatial variability in microtopography and water-table height. 

The goal of this study was to understand the effects of nat¬ 
ural and experimental variation in water-table height on soil 
chemistry, microbial respiratory processes and nutrient cy¬ 
cling in polygonized tundra of the Arctic Coastal Plain in 
Northern Alaska. The motivation was to understand potential 
responses to climate change but also the controls of microto¬ 
pography and water-table height on the current functioning 
of this ecosystem. 

2 Methods 

2.1 Site description and flooding treatment 

This research was carried out in a DTLB in the Barrow 
Environmental Observatory near Barrow, Alaska (71.32° N, 
156.62° W). The DTLB is classified as medium-aged, hav¬ 
ing drained from 50-300 years in the past (Hinkel et al., 
2003). Low-centered polygons have formed over much of the 
basin, particularly in the cental and north sections. Vegeta¬ 
tion is dominated by mosses ( Sphagnum spp.), sedges ( Carex 
aquatilis , Eriophorum spp.), and grasses ( Dupontia fisheri, 
Arctophila fulva), with mosses dominating topographically 
high areas such as polygon rims and Carex dominating low 
areas such as polygon centers (Zona et al., 2009; Olivas 
et al., 2010). The soils within the basin are classified as 
aquahaplels, orthels and histoturbels, with an organic horizon 
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(~12—15 cm thick) overlaying a silty mineral layer. The ac¬ 
tive layer is typically about 30 cm deep (Shiklomanov et al., 
2010). The water table manipulation experiment at this site 
has been described in detail previously (Zona et al., 2009; 
Lipson et al., 2010; Olivas et al., 2010). Briefly, the roughly 
elliptical basin (about 1.5 km long) was divided into three 
sections by plastic dikes buried in the permafrost. The north 
section was designated as the flooded treatment, the central 
as the draining treatment and the south as the control. Water 
was pumped from the central section into the north section. 
However this produced only modest alterations of the wa¬ 
ter table and so additional water was pumped into the north 
section from nearby Cakeeater Lake. This treatment was per¬ 
formed during the summers of 2008 and 2009. The draining 
treatment was very subtle in its effect on water table rela¬ 
tive to the flooding treatment, and the control (south) and 
drained (central) areas were indistinguishable in the response 
variables measured in this study. Therefore the two non- 
flooded areas were pooled for most purposes in this study. 
Raised boardwalks were installed across each of the three 
sections, creating 300 m long transects. All measurements 
presented here occurred along these three transects. Eleva¬ 
tion was recorded with a differential GPS (Trimble R7 and 
R8) every 2 m along each transect, and interpolated linearly 
as needed. For spatially intensive data (see below), eleva¬ 
tion classes were defined by rounding elevation values to the 
nearest 0.1 m. Most of the data presented here were collected 
during the summer of 2009, the second year of the treatment. 
Additional data were collected before the treatment was ini¬ 
tiated (2006-2007), the first treatment year (2008), and the 
summer after the flooding treatment was halted (2010). 

2.2 Soil electrochemistry 

Measurements of pH and oxidation-reduction potential 
(ORP) were taken at 10 cm depth using a portable meter 
(Thermo-Orion). Electrodes were inserted into the soil after 
a small incision was made with a serrated knife. These mea¬ 
surements were made repeatedly at a set of permanent plots 
chosen (based on visual examination of each site) to repre¬ 
sent high and low areas of the landscape within each tran¬ 
sect (three high and three low areas per transect, 18 total). 
These were the same locations used for soil pore water sam¬ 
pling described below. Additionally, more spatially inten¬ 
sive measurements were made on a subset of sampling days. 
ORP measurements are reported as raw values relative to 
the internal Ag/AgCl reference electrode. Dissolved O2 was 
measured using a microelectrode and portable picoammeter 
(Unisense). The O2 electrode was used to find the approxi¬ 
mate location of the transition zone between oxic and anoxic 
horizons (“oxycline”). Because the electrode is only 5 cm 
long, if the oxycline lay deeper than 5cm, it was recorded 
whether the dissolved O2 concentration at 5 cm depth was 
suboxic (<50 % saturation) or oxic (>50 % saturation), and 
assigned non-parametric values 6 or 7, respectively. 
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2.3 Soil water chemistry 

Soil water was collected from 10 cm long, porous soil mois¬ 
ture samplers (Rhizon, Ben Meadows) inserted into the soil 
at an angle, to capture 0-5 cm soil pore water, or vertically, 
to capture water at 5-15 cm depths. The deeper depths were 
only accessible for a limited number of dates later in the sea¬ 
son. The soil moisture samplers were installed in six loca¬ 
tions on each transect, split between high and low areas as 
described above. Water samples were collected through hy¬ 
podermic needles into evacuated tubes (BD Vacutainers). On 
most dates, two sets of water samples were collected; one 
was acidified with a drop of 1M HC1 to preserve the oxida¬ 
tion state of dissolved Fe species. The other set, reserved for 
dissolved gas analysis, was not acidified and was sealed in 
Bitran bags to reduce gas exchange with the atmosphere. All 
water samples were kept refrigerated until analysis. 

Acidified water samples were analyzed for dissolved Fe 2+ 
and Fe 3+ using 1,10-phenanthroline, as described previously 
(Analytical Methods Committee, 1978; Lipson et al., 2010). 
As an index of dissolved aromatic compounds, absorbance 
at 260 nm (A260) was measured on a standardized volume 
(200 pL) in UV-transparent microtiter plates using a mi¬ 
croplate reader (SpectraMax 190, Molecular Devices). Dis¬ 
solved organic C was measured using an assay based on 
Mn(III) as described earlier (Bartlett and Ross, 1988; Lip¬ 
son et al., 2010). 

Ammonium (NH^) was analyzed using a phenolate- 
hypochlorite method and free amino-N, using a ninhydrin 
method (Lipson et al., 1999). Ninhydrin also reacts with 
NH^ so these values were subtracted from ninhydrin-active 
N to yield amino-N. Phosphate (PO4) was analyzed with a 
paramolybdate-ascorbate method (Kuo, 1996). All these as¬ 
says were scaled down for use in a microplate reader. 

The presence of siderophores (biologically produced or¬ 
ganic molecules capable of binding Fe(III) with high affin¬ 
ity) was assayed using a colorimetric assay based on the 
dye, chrome azurol S (Alexander and Zuberer, 1991). This 
dye forms a blue compound when coordinated with Fe(III), 
and so sequestration of Fe(III) from the dye complex by 
siderophores causes a decrease in the blue color. To remove 
interference from dissolved Fe(III) in the water samples, it 
was necessary to dilute samples ten-fold and increase the 
FeCl3 concentration in the dye reagent to 250 pM. Deferox¬ 
amine mesylate (Sigma-Aldrich) was used as a standard. 

Dissolved CH4 and CO2 were analyzed in non-acidified 
water samples. Flasks (25 mL) with butyl stoppers were 
flushed with N2 gas for 1 min, 1.0 mL soil pore water was in¬ 
jected by syringe, allowed to equilibrate at 22 °C overnight, 
and the headspace was sampled by syringe and analyzed 
by gas chromatography (Model 8610C with methanizer and 
flame ionization detector, SRI Instruments, Torrance, CA). 
Original dissolved gas concentrations were calculated from 
the headspace concentrations and Henry’s law constants for 
the two gases. 
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2.4 Microbial activity and biomass 

On three dates in 2009, soil cores were collected to a depth 
of 15 cm using a serrated knife, sealed in two layers of Bitran 
specimen bags, and frozen. Cores were weighed, and the vol¬ 
ume of the frozen cores wrapped tightly in plastic wrap was 
determined by displacement of water. Subsamples of these 
soils (~25 g wet weight) were placed in sealed mason jars 
with septa installed in the lids, and allowed to thaw under 
N2 gas overnight at 4 °C. The following day the jars were 
again flushed with N2 gas to remove accumulated gases, sy¬ 
ringe samples were taken at 6 and 24 h after this point, and 
analyzed by gas chromatography as described above. Other 
subsamples from these cores were subjected to chloroform 
fumigation-extraction for indices of microbial biomass C and 
N. Organic C and ninhydrin-reactive N in these extracts were 
measured as described above. Subsamples were also used 
in proteolysis measurements: soils (~5 g wet weight) were 
placed in tubes with 20 mL sodium acetate buffer (50 mM, 
pH 5.0) with 1 % bovine serum albumin (BSA) as a pro¬ 
tein source, and 0.4 mL toluene to inhibit amino acid up¬ 
take. Samples and soil-free controls were incubated at 6 °C 
on a shaker for 16h. Samples were removed, remaining BSA 
was precipitated by heating at 95 °C for 10 min, and released 
amino acids were measured using the ninhydrin assay. Mi¬ 
crobial biomass and proteolytic rates were also measured on 
soils collected in July 2008. 

Seasonal changes in solid-phase Fe was determined by ex¬ 
tracting frozen soils in 1 N H2SO4 for 24 h. The resulting 
extracts were diluted and assayed for Fe 2+ and Fe 3+ as de¬ 
scribed above. 

A measure of N availability to plants was obtained by 
burying nylon bags with 10 g mixed bed ion exchange resins 
at a depth of 5 cm. The resin was rinsed with copious 
amounts of distilled water before placing in the field, and 
replicate resin bags stored at 4°C served as controls. Resin 
bags were placed in the soil at the beginning of July and col¬ 
lected in mid-August. Resin bags were extracted in 2M KC1, 
analyzed for NHjJ" and NO^, and values from control bags 
were subtracted. This method was used in 2007 (before the 
flooding treatment) and again in 2009. Resin bags were in¬ 
stalled near the same permanent plots referred to above (three 
high and three low locations in each transect). 

2.5 Statistical analysis 

Specific statistical tests used for each data set are described 
in the results section. These included general linear models 
(GLM), analysis of variance (ANOVA), analysis of covari¬ 
ance (ANCOVA), the Kruskal-Wallis test (for analysis of the 
non-parametric oxycline variable), as well as simple regres¬ 
sion and correlation analysis. Systat (v. 12.0) was used for all 
analyses. P-values less than 0.05 are considered significant, 
while 0.1 > P > 0.05 are referred to as marginally significant. 
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3 Results 

3.1 Flooding treatment 

The experimental basin increases in polygonization from 
south to north. This can be seen as a widening range of eleva¬ 
tions from south to central to north treatment areas (Fig. la). 
The north includes areas with elevations as low as 1.5 m, 
but these were omitted from the analysis so that a compa¬ 
rable range of elevations would exist between the flooded 
and reference areas. At any given time, water-table height 
was a function of microtopography, with lower elevation 
areas experiencing larger water-table heights and a dispro¬ 
portionately large effect of the flooding treatment (Fig. la). 
The effects of elevation (F = 103.8, P < 0.001), treatment 
area (N, C and S) (F = 38.0, P < 0.001) and the interac¬ 
tion (F = 29.0, P < 0.001) were all highly significant for 
this ANCOVA (n = 403, R 2 = 0.736). 

In 2008, the flooding treatment initially raised the aver¬ 
age water-table height by 8.6 cm relative to non-flooded ar¬ 
eas, reaching a maximum differential of 15.3 ± 1.2cm by 
the end of summer (data not shown) (Olivas et al., 2010). 
In 2009, the flooding treatment initially achieved an av¬ 
erage water-table height 4.0 ± 1.2 cm higher than the non- 
flooded treatment areas, reaching a maximum differential of 

22.5 ± 1.4 cm at the end of July (Fig. lb). (Late in the sum¬ 
mer of 2009, water was pumped into the south region to 
create an intermediate water table treatment; however, none 
of the data presented here were affected by this treatment). 
The wetter soils conducted heat more efficiently, leading to 
warmer temperatures in mid-summer (Fig. lc). Averaged 
across all dates, temperatures measured in the flooding treat¬ 
ment were 1.0±0.1°C higher than in the unflooded treat¬ 
ments. 

3.2 Soil electrochemistry 

Soil ORP and pH were both controlled by water-table height 
(Fig. 2a, c), with more reducing and less acidic conditions 
in wetter areas. Accordingly, in ANCOVA, ORP was pos¬ 
itively correlated {R 2 = 0.287, F = 68.2, P < 0.001) and 
pH was negatively correlated (R 2 = 0.388, F = 118.6, P < 
0.001) with elevation, and for any given elevation, pH was 
higher (F = 100.1, P < 0.001) and ORP lower (F = 71.9, 
P < 0.001) in the flooding treatment (Fig. 2b, d). The nega¬ 
tive correlation between pH and ORP was highly significant 
(Table 1). 

The seasonal pattern of ORP followed that of water ta¬ 
ble, becoming more oxic as water levels dropped (Fig. 3a, 
b). Despite this seasonal pattern and the overall negative cor¬ 
relation between pH and ORP, pH also increased during the 
early growing season, particularly in low elevation classes 
(Fig. 3c, d), corresponding to a period where reduction of 
solid phase Fe minerals was observed (Lipson et al., 2010), 
see Sect. 3.4 below. 
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Table 1. Correlations among dissolved species in soil water and electrochemical measurements. The upper diagonal shows correlation 
coefficients (with those significant by Bonferroni test, a = 0.05/21 = 0.0024, in bold). Lower diagonal shows sample size for each correlation. 



Fe3 

A260 

DOC 

logCH 4 

logC0 2 

pH 

ORP 

Fe3 


0.824 

0.702 

0.361 

0.561 

0.49 

-0.378 

A260 

126 


0.79 

0.401 

0.593 

0.538 

-0.471 

DOC 

102 

102 


0.292 

0.545 

0.253 

-0.593 

logCH 4 

73 

73 

51 


0.79 

0.505 

-0.46 

logC0 2 

73 

73 

51 

73 


0.602 

-0.553 

pH 

126 

126 

102 

73 

73 


-0.583 

ORP 

90 

90 

68 

55 

55 

720 



Table 2. Oxycline depth index. Values of 0-5 indicate the depth 
in cm at which dissolved O2 dropped to near-anoxic levels. Values 
of 6 and 7 are non-parametric indicators of suboxic (<50 % O2 sat¬ 
uration) and oxic (>50 % O2 saturation) conditions, respectively, at 
5 cm depth. For the seasonal sampling (taken at fewer locations), 
values are means (and standard errors). 



Flooded 

Unflooded 

7/22-7/23 intensive sampling 


median 

3 

7 

mean 

3.61 

6.63 

se 

0.33 

0.09 

n 

49 

75 

Seasonal sampling 


6/29/2009 

4.17(0.91) 

4.17 (0.75) 

7/15/2009 

3.83 (1.45) 

4.93 (0.62) 

7/24/2009 

5 (1.22) 

7(0) 

7/29/2009 

4.17(1.28) 

6.67 (0.26) 


Conditions in the flooded treatment were generally anoxic 
at 5 cm depth or shallower throughout the measurement pe¬ 
riod (Table 2). In the unflooded areas conditions generally 
became oxic in the upper 5 cm by late July, as the water table 
had dropped below 5 cm depth by this time (Fig. lb). The 
effect of flooding on the oxycline depth was highly signifi¬ 
cant for the spatially-intensive 22-23 July data set (Kruskal- 
Wallis test P < 0.001). In the GLM analysis of these data 
(dfe = 120), elevation (F — 48.0, P < 0.001) and flooding 
(F = 54.8, P < 0.001) were also highly significant. For the 
less spatially-intensive seasonal sampling, the flooding effect 
was significant in a two-way ANOVA with date and flooding 
as categorical factors (P = 0.016). 

Trends in pH and ORP disappeared or reversed in the 
year following the water manipulation (2010). Before the 
flooding treatment was initiated and after it ceased, soil pH 
was slightly lower and ORP slightly higher in the north 
(“flooded”) section of the experimental basin (Table 4). 


3.3 Soil water chemistry 

Soil pore water concentrations (from 0-5 cm depth) of aro¬ 
matic compounds (A260), DOC, Fe(III), NH^, amino-N 
and PO4, all declined exponentially throughout the season 
(Fig. 4). Therefore data were log-transformed to linearize 
this trend, and ANCOVA were run using date as a continuous 
variable and topography and flooding as factors (no signifi¬ 
cant interactions with date were found in these models). In 
the case of PO4, the three-way ANOVA (with date coded cat¬ 
egorically) had a significantly lower corrected Aikake infor¬ 
mation criterion (AIC c ) than the ANCOVA, and so these re¬ 
sults are reported instead of the ANCOVA (Table 3). Signif¬ 
icance of effects did not vary between the two model types, 
with one exception (see results for NH^ below). Because 
the effects of flooding and topography were generally inde¬ 
pendent of date, to simplify presentation for most variables 
only the grand means of each treatment category are shown 
in Fig. 5. In the case of NH^, certain dates appeared to 
have a disproportionate effect, and so the seasonal patterns 
for flooded and unfloooded treatments are shown (Fig. 5h). 

A260, DOC, and Fe(III) followed similar trends (Fig. 5a- 
c) and were highly correlated (Table 1). All were signifi¬ 
cantly higher in low elevation areas (Table 3). A260 was sig¬ 
nificantly higher in the flooding treatment; Fe(III) was signif¬ 
icantly higher in low elevation areas of the flooded treatment 
(topography x flooding interaction). The effect of flooding 
on DOC was significant (P = 0.02), and while the response 
appeared to be driven by low elevation areas (Fig. 5b), the to¬ 
pography x flooding interaction was marginal (P = 0.098). 
The trends observed in A260 and Fe(III) during the flood¬ 
ing treatment were reversed in the years before and after the 
flooding treatment was in effect, with generally higher values 
in the areas used for the unflooded treatment (Table 4). 

For two dates later in the season, soil water samples were 
also collected from deeper in the profile (5-15 cm). On these 
dates, dissolved Fe(III) concentrations and A260 values were 
significantly higher at depth, by factors of 2.65 and 2.36, re¬ 
spectively (data not shown). Mean DOC concentrations were 
only 1.33 times higher at depth, and this difference was not 
statistically significant (data not shown). 
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Table 3. Statistical tests (F-statistics and p-values) on the log-transformed water chemistry data (0-5 cm depth) presented in Figs. 5 and 6. 
For all but two variables, ANCOVA models were used with date as a continuous variable (see Fig. 4) and flooding treatment and topography 
as categorical effects. For PO4 and CFI4, three-way ANOVA were used, with date coded categorically (these models had significantly lower 
AIC ( than the ANCOVA). No interactions with date were significant (P < 0.05) in any model. Topo = topographic effect, dfe = degrees of 
freedom in error. Significant effects are shown in bold, marginally significant effects (0.05 < P < 0.1) are italicized. 




Date 

Topo 

Flooding 

Topo x Flood 

dfe 

Fe(Ifl) 

F 

26.127 

92.994 

0.243 

8.345 

121 


P 

<0.001 

<0.001 

0.623 

0.005 


A260 

F 

138.776 

43.739 

18.416 

4.465 

118 


P 

<0.001 

<0.001 

<0.001 

0.037 


DOC 

F 

92.922 

15.874 

5.475 

2.79 

108 


P 

<0.001 

<0.001 

0.021 

0.098 


nh 4 

F 

28.781 

0.03 

1.827 

0.002 

120 


P 

<0.001 

0.864 

0.179 

0.96 


Amino-N 

F 

46.849 

1.781 

0.005 

1.009 

121 


P 

<0.001 

0.185 

0.943 

0.317 


PO4 

F 

248.906 

1.637 

5.551 

1.296 

120 


P 

<0.001 

0.203 

0.02 

0.257 


C0 2 

F 

13.634 

80.015 

0.265 

5.833 

66 


P 

<0.001 

<0.001 

0.609 

0.019 


ch 4 

F 

26.566 

98.053 

0.19 

3.046 

55 


P 

<0.001 

<0.001 

0.664 

0.087 



Table 4. Grand means (and standard errors) in flooding and unflooded treatment areas before (2006-2007), during (2008-2009, shaded) and 
after (2010) flooding treatment. 



pH 

ORP (mV) 

Fe 3 + (pM) 

A260nm (AU) 

year 

unflooded 

flooded 

unflooded 

flooded 

unflooded 

flooded 

unflooded 

flooded 

2006 

5.14 (0.04) 

4.89 (0.07) 

-152.7 (66.2) 

-142.0 (79.3) 

178.6 (8.8) 

165.3 (6.0) 

0.276 (0.019) 

0.238 (0.023) 

2007 

3.86 (0.22) 

3.54 (0.22) 

nd 

nd 

476.1 (15.3) 

337.9(19.9) 

0.308 (0.067) 

0.303(0.129) 

2008 

4.82 (0.08) 

5.00 (0.15) 

180.0 (20.8) 

112.5 (28.5) 

58.9 (7.2) 

74.7 (10.4) 

0.958 (0.217) 

1.380 (0.592) 

2009 

4.62 (0.02) 

5.23 (0.04) 

237.9 (6.4) 

124.3 (10.7) 

34.5 (6.9) 

44.3 (9.6) 

0.224 (0.022) 

0.446 (0.070) 

2010 

5.20 (0.20) 

4.67 (0.22) 

124.8 (34.8) 

220.9 (18.2) 

183.7 (10.6) 

74.9 (8.6) 

0.394 (0.071) 

0.232 (0.071) 


Soil pore water samples were assayed for the presence of 
siderophores on a single day (8 July) in 2009 (Fig. 5d). In 
a two-way ANOVA, the concentration of siderophores was 
much higher in low elevation areas (F = 38.9, P < 0.001) 
and there was no significant effect of flooding. The same 
pattern held in 2010 (data not shown), and including data for 
both 2009 and 2010, siderophore concentrations were highly 
positively correlated with total dissolved Fe (correlation co¬ 
efficient r = 0.931, n = 44) and A260 (r = 0.862, n = 33), 
and negatively correlated with ORP (r = —0.751, n = 26). 

PO4 was significantly higher in the flooded treatment but 
the effect of topography was not significant (Fig. 5e, Table 3). 
There was no significant effect of topography or flooding on 
amino-N in soil pore water (Fig. 5f, Table 3). Neither were 
there consistent effects of flooding or topography on NH^ 


(Table 3). The higher grand means for NHjJ" concentrations 
in the flooding treatment (Fig. 5g) were due to higher val¬ 
ues early in the summer (Fig. 5h), leading to a significant 
date x flooding interaction in the three-way ANOVA with 
date coded categorically (P = 0.031). NOj~ concentrations 
were below detection limits (~5 pM) in the majority of water 
samples, and the mean NO^~ concentrations for flooded and 
unflooded treatments were nearly identical (4.9 vs. 5.1 pM, 
respectively; data not shown). 

Dissolved CH4 and CO2 concentrations were correlated 
with each other and with A260, DOC and Fe(III) (Table 1). 
Like most of the other soil water variables, dissolved CO2 
also declined throughout the season (Fig. 6a). Dissolved 
CH4 followed a different seasonal pattern, peaking in mid¬ 
summer (Fig. 6b). Because of this, the three-way ANOVA 
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6/15 6/29 7/13 7/27 8/10 



Fig. 1. (A) The relationship between water table depth and ele¬ 
vation, and the effects of flooding the North transect (N). Data are 
shown for 21-29 July 2009. Points with elevation below 1.7 m were 
omitted, as this elevation class existed only in the North transect. 
(B) Seasonal changes in water-table height in the three treatment ar¬ 
eas. (C) Soil temperature at 10 cm (measurements associated with 
pH and ORP measurements across each transect, generally made 
between 10 a.m. and 6p.m). 


of dissolved CH4 with date coded categorically had signifi¬ 
cantly lower AIC e than the ANCOVA with date coded as a 
continuous variable (Table 3). At a depth of 0-5 cm, both 
COo and CH4 were found in significantly higher concentra¬ 
tions in topographically low areas (Fig. 6a, b). Across all 
dates, the flooding treatment significantly lowered dissolved 
COo concentrations in high elevation areas (Topo x Flood 
interaction in Table 3, Fig. 6 ); the similar trend for CH4 



Fig. 2. Relationships of oxidation-reduction potential (ORP) and 
pH with water-table height and elevation. (A) ORP vs. water table; 
(B) ORP vs. elevation; (C) pH vs. water table; (D) pH vs. elevation. 
Data are shown for 17-29 July 2009. For the plots showing rela¬ 
tionships with elevation, data were omitted for elevations <1.7m, 
which only occurred in the flooded transect. ORP values are relative 
to a Ag/AgCl reference electrode. 


(Fig. 6b) was marginally significant (P = 0.087, Table 3). 

On one day towards the end of summer (1 August 2009), 
higher dissolved CO2 and CFLj concentrations were observed 
in low elevation areas of the flooding treatment. This led 
to marginally significant date x flooding interactions (P — 
0.067 and 0.063 for CO2 and CH4). This change in the effect 
of the flooding treatment coincided with a period when wa¬ 
ter tables in the control areas were rapidly dropping to near¬ 
minimum levels (Fig. lh); the flooding treatment appears to 
have extended conditions favorable for CH4 and CO2 accu¬ 
mulation later into the season. On this day only (by hap¬ 
penstance), water samples were also collected from 5-15 cm 
depth at the same locations and analyzed for dissolved CO2 
and CH4 (Fig. 6c, d). A three-way ANOVA was run using to¬ 
pography, depth and flooding as factors (dfe = 28). For CH4, 
concentrations were significantly higher at depth {F = 16.1, 
P < 0.001), in low elevation areas (F = 22.9, P < 0.001) 
and in the flooding treatment (F = 11.1, P =0.002), and 
the topography x flooding interaction was also significant 
(F = 5.3, P = 0.029), indicating that the flooding response 
was driven by low elevation areas. For CO2, the effects 
of flooding (F = 8.2, P = 0.008), topography (F = 31.7, 
P < 0.001) and depth (F = 32.7, P < 0.001) were signifi¬ 
cant, as was the topography x flooding interaction (F = 9.4, 
P = 0.005), indicating the differential effect of flooding on 
high vs. low areas, and the depth x topography interaction 
(F = 5.5, P = 0.027), stemming from the steeper depth gra¬ 
dient in topographically high areas. (A steeper gradient along 
the profile implies more rapid diffusion in high elevation ar¬ 
eas). 
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Fig. 3. Seasonal trends in 2009 for (A) ORP by flooding treatment, (B) ORP by elevation class, (C) pH by flooding treatment and (D) pH 
by elevation class. Elevation classes were defined by rounding elevation values to the nearest 0.1 m. All electrochemical and temperature 
readings were taken at 10 cm depth. ORP values are relative to an internal Ag/AgCl reference electrode. 


Table 5. Geometric means (and geometric standard errors) of the CH4:C02 ratio dissolved in soil pore water collected at two depths and in 
emissions from anaerobic laboratory incubations of soils. Means are separated by flooding treatment, microtopography (high vs. low) and 
sampling date. The standard errors are multiplicative, meaning the mean multiplied and divided by the error corresponds to the mean ± 
standard error of the log-transformed data. 


Unflooded Flooded 

Date High Low High Low 


Dissolved 
0-5 cm 


5-15 cm 
Incubations 


6/29/2009 

7/19/2009 

8/1/2009 

8/11/2009 

8/1/2009 


6/29/2009 

7/14/2009 

8/3/2009 


325.4 (1.39) 
90.0 (1.20) 
652.2(1.18) 
485.9(1.19) 
419.2 (1.60) 


957.6 (1.77) 
2118.3(1.29) 
1753.0(1.67) 


136.2(1.33) 
18.1(1.19) 
407.6 (1.56) 
305.6(1.08) 
215.7(1.37) 


360.6(1.58) 
232.5 (1.61) 
372.8 (1.35) 


217.6 (1.09) 

96.3 (1.26) 

508.3 (1.13) 
536.8 (1.15) 
218.0(1.57) 


732.0(1.62) 

713.6 (2.47) 

686.6 (1.25) 


99.4(1.36) 
39.4(1.18) 
116.2(1.05) 
276.7 (1.64) 
59.6(1.50) 


99.5 (1.72) 
30.4(1.29) 
44.7 (1.49) 


Ratios of dissolved CC>2:CH4 were quite high for a wet 
and low oxygen landscape, but generally followed pre¬ 
dictable patterns, decreasing in the most anoxic areas (Ta¬ 
ble 5). This variable was log-normally distributed and so 
values are expressed as geometric means, and a three-way 
ANOVA was run on log-transformed data (in 0-5 cm sam¬ 
ples), with date, flooding and topography as factors. The 


CC>2:CH4 ratio varied by date, dropping in mid-July when 
CH4 concentrations peaked (F = 41.1, P < 0.001), and was 
greater in high areas of the landscape (F = 37.7, P < 0.001). 
Overall, the flooding treatment was not significant, but the 
flooding x date interaction was (F = 2.9, P = 0.042), in¬ 
dicating an inconsistent effect of flooding throughout the 
season. Analyzing each date separately, the increase in 
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Fig. 4. Seasonal patterns of dissolved ion concentration (pM) and 
degree of UV absorbance (absorbance units, AU) in soil pore water 
collected from 0-5 cm depth in 2009. Values are means and stan¬ 
dard errors of both treatments and topographic types. The Y-axis is 
shown on a logarithmic scale. 


CC>2:CH4 ratio in flooded areas on 7/19/2009 and the de¬ 
crease in flooded areas on 8/1/2009 were both marginally 
significant (P = 0.061 and P — 0.05, respectively). On 
8/1/2009 (the day on which both shallow and deeper sam¬ 
ples were analyzed), depth (F = 53, P = 0.028), topography 
(F = 10.1, P = 0.003) and flooding (F = 9.388, P = 0.004) 
were all significant. 

Many soil chemistry variables responded in concert to nat¬ 
ural and experimental changes in water-table height, as seen 
in the intercorrelations of A260, Fe(III), CO2, and CH4; then- 
positive correlations with pH and their negative correlations 
with ORP (Table 1). Furthermore, the logarithm of dis¬ 
solved Fe(III) was positively correlated with the height of 
the water table (r = 0.369, n = 119, P < 0.001), in a re¬ 
gression that incorporated both seasonal and spatial varia¬ 
tion. Interestingly, while Fe(III) and CH4 pore water con¬ 
centrations were positively correlated (because both increase 
with anoxia), the logarithm of dissolved CH4 concentration 
was negatively correlated with the ratio of Fe 3+ /total Fe 
(Fig. 7; R 2 = 0.406, F = 14.2, P < 0.001). This relation¬ 
ship was significant while controlling for the effect of topog¬ 
raphy (F =21.8, P < 0.001), and the slopes were not sig¬ 
nificantly different between high and low topographic types 
(F = 1.598, P = 0.21). This indicates that this negative re¬ 
lationship was not mediated solely by O2 (inasmuch as the 
topographic types control for variation in O2), and so a di¬ 
rect relationship might exist between dissolved CH4 and the 
average oxidation state of the soluble Fe pool. 

3.4 Microbial activity and biomass 

Consistent with the effects of flooding and microtopogra¬ 
phy on ORP and DO shown above, anaerobic microbial pro¬ 
cesses also responded to natural and experimental variations 
in water-table height (Fig. 8). Potential CO2 and CH4 pro- 


Flooded Unflooded Flooded Unflooded 



Fig. 5. Soil pore water chemistry extracted from upper 5 cm of soil 
in summer, 2009. (A-G) Grand means of all samples dates in sum¬ 
mer 2009 in topographically high and low areas of flooded and un¬ 
flooded treatments. Siderophore concentrations are presented as pM 
deferoxamine equivalents. (FI) Seasonal trends for NH^ in flooded 
and unflooded treatments. Refer to Table 3 and the text for statisti¬ 
cal significance. 


duction rates were measured in anaerobic laboratory incuba¬ 
tions on three dates, and effects were tested using a three- 
way ANOVA with date, topography and flooding as factors. 
CO2 production rates in anaerobic incubations were high¬ 
est on the first date (P = 0.003), and consistently greater in 
high topography regions of the flooding treatment (Fig. 8a, 
topography x flooding interaction: P = 0.019). The direct 
topographic (P = 0.018) and flooding (P = 0.003) effects 
were also significant, but both appear to be driven by high 
areas in the flooding treatment. In other words, the poten¬ 
tial for anaerobic respiration mainly responded to flooding 
in the high regions, presumably because microbial commu¬ 
nities in low regions were already predisposed to anaerobic 
respiration. CH4 production rates in these incubations were 
greater in low elevation areas and in the flooding treatment 
(P < 0.001 for both; date and interactions not significant. 
Fig. 8b). The C02:CH4 ratios evolved in incubations were 
similar to those dissolved in soil pore water, and showed 
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Fig. 6. Dissolved gases in soil pore water from high and low areas 
within flooded and unflooded treatments: (A) CO2 and (B) CH4 at 
a depth of 0-5 cm over four dates in 2009; (C) C0 2 and (D) CH 4 
at 0-5 cm and 5-15 cm depth on 1 August 2009. The Y-axes are 
scaled logarithmically. Refer to Table 3 and the text for statistical 
significance. 



Fig. 7. Relationship between concentration of CH4 (log- 
transformed) and the ratio of Fe 3+ to total Fe dissolved in soil pore 
water in topographically high and low areas. Separate regression 
lines are shown for low (solid) and high (dashed). 


similar spatial patterns (Table 5). C0 2 :CH 4 was reduced in 
flooded areas (F = 19.5, P < 0.001) and in topographically 
low areas (F = 51.8, P < 0.001). Differences among dates 
were not significant in this analysis. 

Dissimilatory Fe(III) reduction over the course of the sum¬ 
mer was monitored as changes in the reduction state of Fe- 
bearing minerals in the soil. Acid-extractable Fe(III) de¬ 
clined throughout the season as Fe(II) increased, indicat¬ 
ing net Fe reduction (Fig. 8c). There was significantly 
more acid-extractable Fe(II) in low areas of the landscape 
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Fig. 8. Indicators of anaerobic processes in soils collected in 2009. 
Production rates of (A) C0 2 and (B) CH4 in anaerobic laboratory 
incubations of soil from flooded and non-flooded treatments. (C) 
Seasonal changes in acid-extractable Fe(III) and Fe(II) from high 
and low areas of the landscape. 


(P = 0.015). Flooding effects were not statistically signifi¬ 
cant. Seasonal Fe reduction rates calculated from either the 
increase in Fe(II) or the decrease in Fe(III) in the acid solu¬ 
ble fraction were greater in low elevation areas (high: 6.4- 
27.4 mmole Fe m -2 d -1 ; low: 28.0-36.4 mmole Fem _2 d _1 , 
with low and high ranges based on AFe(II) and AFe(III), re¬ 
spectively). 

In contrast to these measures of anaerobic activity, more 
general measures of microbial biomass showed more sub¬ 
tle responses to water table. Microbial biomass C, mea¬ 
sured by fumigation-extraction, showed no significant ef¬ 
fects of flooding, topography or date (Fig. 9a). Microbial 
biomass N (measured as ninhydrin-active N released by fu¬ 
migation extraction) showed marginally significant effects 
driven by lower means in flooded, low-elevation areas in 
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Fig. 9. Measures of microbial biomass and N cycling in soils col¬ 
lected from high and low areas of flooded and unflooded areas: (A) 
Microbial biomass C and (B) ninhydrin-active N measured by the 
fumigation-extraction method, (C) potential proteolysis rates, and 
(D) NH^" bound to buried resin bags over the 2007 (pre-flooding) 
and 2009 (post-flooding) seasons. 


2009 (Fig. 9b). In a GLM analysis (dfe = 44), date was 
significant (F = 9.106, P < 0.001), the flooding effect (F = 
3.252, P = 0.078) and the topography x flooding interac¬ 
tion (F = 3.215, P — 0.080) were marginally significant, and 
the direct topography effect was not significant (F = 1.687, 
P= 0.201). 

Soil protease activity varied by date (F = 0.016), but 
flooding and topographic effects were not significant 
(Fig. 9c). In contrast, N availability integrated over the sum¬ 
mer appeared to respond to water-table height: more N was 
captured by buried resin bags in low elevation areas dur¬ 
ing the pre-flooding year (2007) and in the flooding treat¬ 
ment in 2009 (Fig. 9d). In the three-way ANOVA comparing 
year, treatment area, and topography (dfe = 30), topography 
(F = 4.421, P = 0.044) and the year x transect interaction 
(F = 5.218, P = 0.030) were significant. (This interaction 
indicates the effect of flooding, which occurred in 2009, but 
not 2007). 


4 Discussion 

Being a relatively flat area with drainage blocked by per¬ 
mafrost and a climate with generally low evaporative de¬ 
mand, this ecosystem is subject to regular flooding from 
snow melt. Therefore, experimental flooding mainly accen¬ 
tuated the normally occurring responses to seasonal and spa¬ 
tial variation in water-table height. Microtopography aris¬ 
ing from ice wedge polygons controlled spatial variation in 
water-table height, leading to sharp contrasts in soil chem¬ 
istry and microbiology over short spatial scales. The ef¬ 
fects of the flooding treatment recovered quickly after the 
water table manipulation ended: during the flooding treat- 


Seasonal Flooding Micro¬ 
variation \ I ✓ topography 

+ Water Table 


- Gas diffusion rates 


/ 

- Dissolved 0 2 

* \ 

+ ligands 

( A26 . 0nm > / (production/ 

I concentration) 

+ dissolved Fe(lll) 

l 

+ pH 


Fig. 10. Synthesis of major results from this paper, illustrating pos¬ 
sible mechanisms by which water table affected the measured vari¬ 
ables. The dotted flat arrow from Fe(III) to CH4 indicates a possible 
negative effect. 


ment ORP was lower and pH, Fe(III) and A260 were higher 
in the flooded section compared to the unflooded sections, 
but these trends reversed in 2010 after the flooding treatment 
had ended. In the absence of the experimental manipulation, 
the less polygonized southern end of the basin was more re¬ 
ducing, on average, than the more polygonized northern sec¬ 
tion. Increased polygonization leads to a greater abundance 
of elevated, oxic areas within the landscape, which raised the 
overall average ORP of the landscape. 

Many of the observed effects of natural and experimental 
variation in water-table height can be attributed to lower dif¬ 
fusion rates of gases in wet soils (Fig. 10). Inundated areas 
become anoxic, and the microbial communities appear to re¬ 
spond by producing organic ligands that solubilize Fe(III). 
The process of solubilizing and reducing Fe(III) (as well as 
other alternative e _ acceptors) consumes protons, leading to 
higher pH (Reddy and DeLaune, 2008). Conversely, protons 
are released and pH is lowered under oxidizing conditions; 
this was seen dramatically in 2007 (Table 4), an exception¬ 
ally dry year in which water-table height dropped to much 
lower levels then in other years of the study (Olivas et al., 
2010). Flooded, anoxic conditions also promote the growth 
and activity of methanogens, and lower diffusion rates can 
also trap more CH4 in the soil. However, the increased avail¬ 
ability of Fe(III) could competitively inhibit methane pro¬ 
duction (Roden and Wetzel, 1996; Bond and Lovley, 2002; 
Kiisel et al., 2008; Jerman et al., 2009), weakening these 
effects (dotted flat arrow in lower left, Fig. 10). High wa¬ 
ter tables also slow diffusion of COo out of the soil, lead¬ 
ing to higher concentrations in lower, wetter areas. It is also 
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possible that microtopographic patterns of plant growth and 
litter accumulation contribute to higher dissolved CO2 con¬ 
centrations in polygon centers (dotted arrow in upper right, 
Fig. 10) (Zona et al., 2011). Root respiration from vascu¬ 
lar plants such as Carex aquatilis, which are more preva¬ 
lent in the lower areas, could also contribute to this trend. 
Conversely, submerged stems of Sphagnum spp., dominant 
in higher parts of the landscape, may efficiently draw down 
dissolved CO2 (Rydin and Clymo, 1989). 

The inference that microbes solubilize Fe(III) in response 
to anoxia is based on (1) the surprisingly high concentrations 
of Fe(III) in soil pore water, given only moderately acidic 
conditions, best explained by the presence of organic chelat¬ 
ing agents; (2) the high concentrations of siderophores in 
pore water of low, anoxic areas; and (3) the strong correla¬ 
tions among soluble Fe, siderophores, the UV absorbance of 
soil pore water (A260) and ORR Aromatic compounds, in¬ 
cluding several classes of siderophores (Drechsel and Jung, 
1998), strongly absorb UV wavelengths. The presence of Fe 
can also enhance the UV absorbance of solutions (Weishaar 
et al., 2003), though the correlations of A260 with DOC and 
siderophores show that the relationship between Fe(III) and 
A260 was not produced entirely by this effect. 

It has already been established that the process of Fe(III) 
reduction is important in this ecosystem (Lipson et al., 2010). 
However, the present study shows the importance of this pro¬ 
cess in shaping fine scale patterns of soil chemistry. This 
is expressed in the solubilization of Fe(III) in response to 
anoxic conditions, and the inverse microtopographic pat¬ 
terns of redox conditions and pH, indicating that proton¬ 
consuming reactions such as Fe(III) reduction take place 
where reducing conditions predominate. Furthermore, it 
was not known how the prevalence of Fe(III) reduction im¬ 
pacts CH4 production in this ecosystem. The high ratios 
of CC>2:CH4 dissolved in soil pore water and produced in 
anaerobic laboratory incubations indicate an ecosystem with 
an abundance of alternative e _ acceptors for anaerobic res¬ 
piration (Updegraff et al., 1995; van Hulzen et al., 1999; 
Keller and Bridgham, 2007; Keller et al., 2009). In the ab¬ 
sence of alternative e _ acceptors, roughly equimolar con¬ 
centrations of CO2 and CH4 would be produced (i.e., a ra¬ 
tio of one), whereas values observed in incubations and dis¬ 
solved in soil pore water were generally orders of magni¬ 
tude higher. Hence, even in the most anoxic areas of the 
landscape, methanogenesis still constituted a minor fraction 
of metabolic e _ flow (at least for the surface layers stud¬ 
ied here). The negative relationship between dissolved CH4 
concentrations and the oxidation state of the soluble Fe pool 
(Fe 3+ /total) shown in Fig. 7 provides further evidence that 
methanogenesis is competitively inhibited by alternative e - 
acceptors. O2 is probably the primary control over dissolved 
CH4 in the surface layers, as can be seen in the lower CH4 
concentrations and higher CC>2:CH4 ratios in raised, more 
oxic areas of the landscape. Likewise, the average oxida¬ 
tion state of Fe is higher in oxic areas. However, the same 
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negative relationship between CH4 and Fe 3+ /total exists for 
both high and low areas, indicating a direct relationship be¬ 
tween these variables rather than one that is purely mediated 
by O2. The oxidation state of the dissolved Fe pool could be 
interpreted as a measure of Fe(III) availability for dissim- 
ilatory Fe reduction, as the tendency for Fe(III) reduction 
given by the Nemst equation depends on the ratio of oxi¬ 
dized and reduced Fe. However, it is also true that chelated 
Fe reacts rapidly with soil humic substances and so various 
Fe(III) chelates have been used to measure the oxidation state 
and “reduction capacity” of humic acids in laboratory incu¬ 
bations (Peretyazhko and Sposito, 2006). Our soils are novel 
in having high concentrations of chelated Fe(III) in soil pore 
water, and so the oxidation state of the dissolved Fe pool 
might provide an in situ measure of the oxidation state of 
the humic substances in contact with this soluble Fe pool. 
If this is the case, then the negative relationship shown in 
Fig. 7 might be caused by preferential use of oxidized humic 
substances as terminal e“ acceptors for anaerobic respiration 
over methanogenesis, as seen in other studies (Cervantes et 
al., 2000; Heitmann et al., 2007). Addition of the humic ana¬ 
log, AQDS, stimulated respiration in laboratory incubations 
of soils from this ecosystem (Lipson et al., 2010), showing 
that quinones found in humic substances can serve as termi¬ 
nal e _ acceptors in these soils. 

The availability of alternative e“ acceptors (Fe(III) 
and/or humic substances) and the high CC>2:CH4 ratios 
explain why reported CH4 fluxes from this landscape 
(24.6mgCH4 m -2 d -1 ) (Zona et al., 2009) are near the low 
end for those reported in the Arctic (Blodau, 2002; Wagner 
et al., 2003). And so in terms of climate change, the poten¬ 
tial for increased methane production from ecosystems such 
as the DTLB studied here could be limited by the high avail¬ 
ability of Fe(III) and humic substances. This inhibitory effect 
could be enhanced if deeper thaw leads to mobilization of Fe 
from the mineral layers that currently reside in the upper per¬ 
mafrost of many DTLB’s of the Arctic Coastal Plain (Hinkel 
et al., 2003). 

The results of this study indicate that microbial commu¬ 
nities in much of this landscape are dominated by faculta¬ 
tive and/or strict anaerobes. It was mainly in high elevation 
areas (e.g., polygon rims), where conditions were normally 
oxic near the surface for most of the year, where the po¬ 
tential for anaerobic respiration responded to the flooding 
treatment (Fig. 8a), indicating that anaerobically-respiring 
microbes were already abundant in low areas of the land¬ 
scape. Methanogens, being strict anaerobes that thrive in the 
most reducing conditions, responded more dramatically to 
the flooding treatment and showed a marked preference for 
low elevation areas, where the risk of exposure to O2 is lower 
(Fig. 8b). This landscape pattern has been noted in other 
studies (Wagner et al., 2003; Strack and Waddington, 2007). 
There were only the subtlest of hints that flooding might cre¬ 
ate suboptimal conditions for microbial activity: lower dis¬ 
solved CO2 concentrations in flooded, topographically high 
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areas could have been caused by a depression of respiration 
due to suboxic conditions (Fig. 6a). Microbial biomass C 
was not impacted by flooding, though one could argue that 
the statistically weak decrease in microbial biomass N in low 
elevation areas of the flooding treatment represented a nega¬ 
tive microbial response to flooding. However, there was no 
consistent trend toward lower microbial biomass C or N in 
lower elevation areas, indicating that anoxic conditions do 
not generally limit the biomass of the local microbial com¬ 
munity. Flooding experiments in northern peat soils have 
often produced insignificant or relatively minor effects on 
ecosystem respiration (Updegraff et al., 2001; Chivers et al., 
2009; Huemmrich et al., 2010), while lowering water table 
tends to have more dramatic effects (Oechel et al., 1998; 
Strack et al., 2006; Stack and Waddington, 2007; Som- 
merkorn, 2008). This dichotomy is understandable given the 
tendency of these ecosystems toward waterlogged, suboxic 
conditions. Increases in soil CO2 fluxes with lowered wa¬ 
ter tables is generally attributed to increased O2 availabil¬ 
ity (Billings et al., 1982; Sommerkorn, 2008; Olivas et al., 
2010). While O2 certainly plays a major role, the impressive 
capacity for anaerobic respiration of some peat soils should 
not be overlooked. At our study site, the flooded treatment 
area was submerged under water for nearly the entire 2008 
growing season, and yet chamber-based measurements of 
ecosystem respiration (ER) were only 23 % less than those 
measured in the control area, where the water table was gen¬ 
erally below the surface (Olivas et al., 2010). Moreover, 
landscape scale measurements of net ecosystem CO2 ex¬ 
change and estimation of ER from eddy covariance showed 
that the flooded area lost CO2 more quickly with higher wind 
speed (Zona et al., 2012). Similarly, eddy covariance mea¬ 
surements of CH4 fluxes from this ecosystem also showed 
strong effects of wind speed (Sturtevant et al., 2011). These 
wind-driven fluxes would not be captured by chamber mea¬ 
surements and so the decrease in ER reported with flooding 
(Olivas et al., 2010) may have been over-estimated. Our ob¬ 
servation in the current study of high dissolved CO2 concen¬ 
trations in low, wet areas and steeper CO2 concentration gra¬ 
dients in higher, drier areas suggests that the effects of soil 
water on diffusion rates of CO2 should be considered when 
interpreting the effects of water table on measurements of 
ecosystem respiration. 

The increased concentrations of P seen in the flooded treat¬ 
ment could be caused by redox-mediated release of P oc¬ 
cluded in ferric hydroxides, as has been observed in other 
studies (Chambers and Odum, 1990; Shenker et al., 2005). 
However, if this were the case, one would expect to see 
higher dissolved P in topographically low areas. The lack 
of response of soluble P to topography makes this explana¬ 
tion less likely. The nutrient cycling effects might also be 
explained by the temperature/thaw depth aspect of the flood¬ 
ing treatment. The largest increases in nutrient concentra¬ 
tions were shortly after break up, before soil temperatures 
had warmed in response to the flooding treatment. The in- 


589 

creases in N and P concentrations observed in the flooding 
treatments early in 2009 could have been partly caused by 
the release of labile substrates from permafrost due to the in¬ 
creased thaw depth in 2008 (Shiklomanov et al., 2010). The 
increased N bound to buried resin bags could have also have 
been a result of increased soil temperatures in the current 
year leading to greater rates of N mineralization. In high Arc¬ 
tic regions of continuous permafrost, where thermokarsting 
is leading to expansion of flooded areas (Smith et al., 2005), 
warming and flooding could have a positive interactive ef¬ 
fect on nutrient availability. This flush of nutrients could 
stimulate plant production, partially offsetting C losses from 
permafrost (Shaver et al., 1992), but it could also ultimately 
stimulate decomposition and further loss of soil C (Mack et 
al., 2004). 

The observed pulse of nutrients after snow melt is com¬ 
monly reported in arctic, alpine and other cold ecosystems 
(Chapin et al., 1978; Zak et al., 1990; Brooks et al., 1998; 
Lipson et al., 1999; Edwards et al., 2006). Maximum thaw 
depths occur in late summer after plants have already begun 
to senesce (Olivas et al., 2010) and so nutrients made avail¬ 
able from thawing of permafrost might not be accessible to 
plants until the following break up. The rapid decline after 
snow melt of dissolved nutrients probably represents rapid 
uptake by plants, as drainage from the basin is blocked by 
permafrost and little water is lost to run-off (Hinkel et al., 
2001). This early season pulse has long been hypothesized to 
represent a major contribution to plant nutrient requirements 
in Arctic coastal tundra (Chapin et al., 1978), and such a pat¬ 
tern would help plants capitalize on nutrients mobilized from 
recently thawed permafrost. 

5 Summary and conclusions 

The flooding treatment produced two distinct suites of re¬ 
sults: (1) elevated water-table heights led to more reducing 
conditions, stimulating siderophore production and the sol¬ 
ubilization and reduction of Fe(III); (2) Flooding stimulated 
N and P availability, probably due to warmer temperatures 
and deeper active layer depths. This study provides evi¬ 
dence that methanogenesis is inhibited in this ecosystem by 
the availability of Fe(III) and humic substances as alterna¬ 
tive e _ acceptors for anaerobic respiration. This study also 
demonstrates how dissimilatory Fe(III) reduction interacts 
with ice wedge polygon-induced microtopography to shape 
landscape patterns of soil chemistry and microbiology. 
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